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Abstract 

In this paper, we formalize the notion of A-AT-model (where A is a non- 
null integer) for a given chain complex, which allows the computation of ho- 
mological information in the integer domain avoiding using the Smith Normal 
Form of the boundary matrices. We present an algorithm for computing such 
a model, obtaining Betti numbers, the prime numbers p involved in the in- 
variant factors of the torsion subgroup of homology, the amount of invariant 
factors that are a power of p and a set of representative cycles of generators 
of homology mod p, for each p. Moreover, we establish the minimum valid 
A for such a construction, what cuts down the computational costs related 
to the torsion subgroup. The tools described here are useful to determine 
topological information of nD structured objects such as simplicial, cubical 
or simploidal complexes and are applicable to extract such an information 
from digital pictures. 

Keywords: algebraic topological model, nD digital image, integer 
homology, chain complex 



1. Introduction 

Important questions in Vision and Image Processing involve an accurate 
topological analysis. In some scientific areas, such as Medical Imaging and 
Remote Sensing, there is a growing need to visualize, analyze and manipu- 
late data of high complexity and dimensionality. We are concerned about 
providing algorithms to solve topological problems in any dimension. For 
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this aim, efficient methods from combinatorial or algebraic topology to com- 
pute topological features are needed. Due to the fact that we will work 
with n-dimensional combinatorial objects (mainly, n = 3,4...), topological 
features are not limited to Betti numbers, Euler characteristic, connectivity, 
number of holes or cavities, but also include other advanced characteristics 
such as (co)homology in the integer domain, cohomology ring, cohomology 
operations or homotopy groups, which can help to discriminate topologi- 
cally non-equivalent objects. Algebraic Topology and Homological Algebra 
fields provide the necessary tools to capture all this topological stuff (see 



.14,1121,11; 

The classical algorithm for computing homology groups (in the integer 
domain) uses the Smith Normal Form (SNF) of the boundary matrices of a 
chain complex. Explicit examples can be given for which the computation of 
SNF has a worst-case computational complexity which grows exponentially 
in both space and time l3|. Many algorithms have been devised to improve 
this complexity bound [ill, 18, lil 15 



In the AT- model theoryp, |6[ , starting from a simplicial complex K, a 
chain homotopy equivalence is generated from the chain complex associated 
to K, C{K), to its homology. In this case, the complexity of the associated 
algorithm is 0{m^), where m is the number of generators of the complex, 
but the computation is carried out over a field. In this paper, we extend 
the notion of AT-model for integer homology computation. More concretely, 
we define the homological algebra notion of A-AT-model and describe an 
algorithm for computing a structure of this kind from any chain complex, 
with coefficients in the integer domain (Z). Starting from this information, 
it is possible, in particular, to obtain Betti numbers, the prime numbers p 
involved in the invariant factors (corresponding to the torsion subgroup of 
the homology), the amount of invariant factors that are a power of p and a 
set of representative cycles of generators of homology mod p, without using 
SNF. We present an algorithm for computing such a A-AT-model for a given 
chain complex and extracting all this homological information, performed in 
0{m^ip{X)) in the worst case, ip being the Euler function. 

In the following section, we recall classical definitions from Algebraic 
Topology. In Section 3, previous tools for computing topological information 
such as AM-models (which make use of SNF), and AT-models (for computing 
homological information over a field) are given. In Section 4, we define the 
notion of A-AT-model, study its properties, give an algorithm for computing 
it and study its complexity. Finally, we describe how to extract homological 
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information in the integer domain from A-AT-models. The last section is 
devoted to conclusions and future work. 



2. Background 

This section introduces the basic concepts and definitions needed through- 
out the paper, taking as main reference Munkres' book 14 



A chain complex C is a sequence ■ ■ ■ ^ Cq Cq-i ^ ■ ■ ■ 
of free abelian groups Cq and homomorphisms dq : Cq ^ Cq-i, such that, for 
all q, dqdq+i = 0. An element a G Cg is a q- chain of C; we say that q is 
the dimension of a, and write q = dim{a). The set d = {dq} is called the 
differential of C. If there is no confusion, we can write a G C if a G Cg for 
some q. Similarly, we can omit the subindex q in dq{a). 

The chain complex C is finite if there exists an integer n > such that 
Cg = for g > 77, and each abelian group Cq is finitely generated. In this case, 
if C„ 7^ 0, we say that dim of C is n. 

Since our goal is the computation of homological information of finite 
objects, all chain complexes considered here are finite. In this case, C can be 
encoded as a pair (C, d), where: (1) C = {Cq}, and Cq is a set of generators 
of Cq] (2) d = {dq} and for each q, dq is the differential of C in dim q with 
respect to the bases Cq and Cg_i. Sometimes, we write C = {C,d) if the 
chain complex C is encoded as the pair {C,d). If b ^ Cq and a G Cq, then 
{b, a) denotes the coefficient of the generator a in the g-chain b when b is 
expressed linear combination of elements of Cq. 

Example 2.1. Shapes are classically modeled with cellular subdivisions. Sev- 
eral combinatorial structures may represent such a subdivision. Simplicial 
complexes have proved to be a useful tool to model a geometric object. Roughly 
speaking, they are collections of simplices ( vertices, edges, triangles, tetrahe- 
dra, ...) that fit together in a natural way to form the object. A q-simplex 
a is a convex hull of a set of q affinely independent points. If the set is 
{vq, vi, . . . , Vq} then the simplex is denoted by VqVi ■ ■ - Vq. The boundary oper- 
ator over the q-simplex a is defined by dqiv^vi ■ ■ - Vq) = ^^=o(~-^)*'^o'yi ■ ■ - Vi- ■ - Vq 
where Vi means that the vertex Vi is omitted. For every simplicial complex 
K, one can define a chain complex C{K) (canonically associated to it) such 
that the q-chains are linear combinations of q-simplices. By linearity, the 
boundary operator can be extended to chains to define the differential d of 
C{K). The homology of K is then defined as the homology ofC{K). Another 
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Figure 1: The Torus and a triangulation T of it. 
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Figure 2: The Klein bottle and a triangulation of it. 



way of combinatorial representation of a geometric structure (which could 
arise naturally, for example, from tomography, numerical computations and 
graphics ), is by means of cubical grids, which subdivide the space into cubes 
with vertices in an integer lattice. This approach, that can be generalized to 
an arbitrary dimension, leads to a cubical complex. The homology of a given 
cubical complex is the homology of the cubical chain complex associated to 
it Finally, simploidal sets 0/ include simplicial complexes and cubical 
complexes as particular cases. They can be used for representing 'hybrid' 
grids coming from finite element methods. In lid], a free chain complex is 
associated to a simploidal set and the homology of the simploidal set is defined 
as the homology of the associated chain complex. 



Given a chain complex C with differential d, a g-chain a G Cg is called 
a q-cycle if dq{a) = 0. If a = dq^i{b) for some b G Cg+i then a is called a 
q-boundary. Denote the groups of g-cycles and g-boundaries by Zg and Bg 
respectively. We say that two g-cycles a and b are homologous if there exists 



a (g + l)-chain c such that a = b + dg^i{c). Define the integer qth homology 
group to be the quotient group Zg/Bg, denoted by Hg{C; Z). 

For each g, the integer gth homology group Hg{C] Z) is a finitely generated 
abelian group. Moreover, Hg{C; Z) is isomorphic to Fg{C; Z) ©Tg(C; Z) where 

FgiC; Z) = Z © ■ ■ ■ © Z and T,(C; Z) = (Z/a(,,i)) © ■ ■ ■ © (Z/a(,,,)) 

are the free subgroup and the torsion subgroup of Hg{C; Z), respectively. The 
rank of Z), denoted by is called the gth Betti number of C. In- 
tuitively, /3o is the number of connected components, f3i is the number of 
independent holes and is the number of cavities. Each is a power 

of a prime, = P*^''^'^ They are called the invariant factors of Hg{C; Z). 
The numbers /3q and are uniquely determined by Hg{C; Z) (up to a 

rearrangement). Therefore, this representation is in some sense a 'canoni- 
cal form' for Hg{C; Z). For all q, there exists a finite number of elements 
of Hg{C; Z) from which we can deduce all the Hg{C; Z) elements. Those 
elements are called homology generators of dim q. We say that a is a repre- 
sentative q-cycle of a homology generator a of the free subgroup of Hg{C; Z) 
if a = a + Bg and for each fi E Z, fi ^ 0, fia ^ Im d. We say that a is a 
representative q-cycle of a homology generator a of the torsion subgroup of 
Hg{C] Z) if a = a + Bg and there exists yU G Z, /i 7^ 0, such that fia G /md. 
We denote a = [a]. 

Proposition 2.2. fl^, p. 332] The qth homology group ofC with coefficients 
in Z/p forp being a prime, denoted by Hg{C; Z/p), is a vector space. Its rank, 
denoted by P(q,p), depends on the prime p and it is related to the number of 
invariant factors of Hg{C; Z) that are a power of p, T(^q,p), via the Universal 
Coefficient Theorem for Homology, which implies that, for each prime p, 

T(g,p) = /3(g,p) - Pg- T(^q~i,p) for q > and T(o,p) = /3(o,p) - /3o- 

Example 2.3. Consider a triangulation T of the torus and a triangulation 
S of the Klein bottle (see Fig. [ip. Let C{T) and C{S) be the chain complexes 
canonically associated to T and S, respectively. 
The homology groups of the torus are: 

HoiCiT); Z) ^ Z, i/i(C(T); Z) ~ Z © Z, H^iCiT); Z) ~ Z. 

The Betti numbers for the torus are: /3q = 1, /3f = 2 and /3J = 1. Repre- 
sentative cycles of homology generators of dim 1 are shown in red and green 
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in Fig. respectively. The 2-chain r which is the sum (up to signs) of all 
the triangles ofT satisfies that 92(r) = 0, so t is a representative cycle of a 
homology generator of dim 2. 

In the case of the Klein bottle, its homology groups are: 

Ho{C{S); Z) ~ Z, H^{C{S); Z) ~ Z © Z/2, H^iCiS); Z) ~ 0. 

The Betti numbers of the Klein bottle are = 1, PI = 1 and /3| = 0. A rep- 
resentative cycle of a homology generator of the free subgroup of Hi{C{S); Z) 
is shown in red in Figl^ The invariant factor of Hi{C{S);Z) is 2. The 
representative cycle a of a homology generator of the torsion subgroup of 
Hi{C{S); Z), which is shown in green in Figure\^ satisfies that 92(7) = 2a 
where 7 is the sum (up to signs) of all the triangles of S. 

Working with coefficients in a field, it is not possible, in general, to get 
all the information corresponding to the torsion subgroup of the homology. 
For example, the homology groups of S with coefficients m Z/3 are 

HoiCiS); Z/3) = Z/3, H^{C{S); Z/3) = Z/3, H^iCiS); Z/3) = 

since (92(7) = 2a, so ^2(27) = a. Notice that the torsion part is lost since 
the prime p = 3 is not involved in any invariant factor of H{C{S), Z). The 
homology groups of the Klein bottle with coefficients in Q are 

HoiCiS); Q) = Q, //i(C(5); Q) = Q, H^iCiS); Q) = 

since 52(^7) = a. Finally, the homology groups of the Klein bottle with 
coefficients in Z/2 are 

Ho{C{Sy, z/2) = z/2, H,{C{Sy, z/2) = z/2 © z/2, H2{C{Sy, z/2) = z/2 

since ^2(7) = 2a = 0. 

Let C and C be two chain complexes. Denote the differential of C and C by 
d and d', respectively. A chain map / : C — )■ C is a family of homomorphisms 
{fq '■ Cq — )■ Cq} such that d'^fg = fq-idq for all g > 0. If no confusion can 
arise, we omit the subindex and write /(a) instead of fq{a) for any q and 
a E Cq. A chain map f : C ^ C induces a homomorphism [/] : H{C] Z) — >■ 
H{C'-Z) where [/]([a]) = [/(a)] for every [a] e H{C;Z). U f,g : C ^ C 
are chain maps, then a chain homotopy : C — ?■ C of / to is a family of 
homomorphisms {(f)q : Cq — )■ C'^^^} such that fq — Qq = d'^j^^^cpq + (pq^idq. A 
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Figure 3: The chain contraction {f,g,4') of C[K) to C{L). 



chain contraction of a chain complex C, to another chain complex C with 
differentials d and d! , respectively, is a set {f,g,4>) such that: f : C ^ C 
and g : C ^ C are chain maps; fg is the identity map of C (denoted by 
idc') and : C — C is a chain homotopy of the identity map of C to gf, 
that is, (f)d + d(f) = idc — gf - Important properties of chain contractions are: 
(1) g is injective and / is surjective since fg = idc' and, in particular, C 
has fewer or the same number of generators than C; (2) [f][g] = idn^c'-.z) 
and [g][f] = idH(C:z) and, therefore, the homology groups of C and C are 
isomorphic (l3| . 

Example 2.4. Let K he the simplicial complex given by the following set of 
simplices (see Fig. on the left): 

{Vl, V2, f3, V4, ViV2, V2V3, V2V4, V1V2V4, V2V3V4} 

and let L be the simplicial complex given by the vertex v^. 

Let f : C{K) — > C{L) and g : C{L) — )■ C{K) be the chain maps given by 
f{vi) = V4 for 2 = 1, 2, 3, 4 and f is null over the other simplices; g{v4) = V4. 
Let (j) : C{K) — )■ C{K) be the chain homotopy ofidc^K) to gf given by: 0(f 1) = 
V1V4, 0(^2) = V2V4, (t>{v3) = V3V4, (t>{viV2) = V1V2V4, (f){v2V3) = V2V3V4 and 4> is 
null over the other simplices. The set (/, g, (p) is a chain contraction ofC{K) 
to C{L) and hence, C{K) and C{L) have isomorphic homology groups. Fig. 
illustrates such a chain contraction. Observe that the action of the map (j) 
'reduces ' the simplicial complex to its homology. 
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3. AM-model and AT-model 

We recall here the concepts of AM-model and AT-model. Both tools 
represent the prior work that lays the foundations of A- AT-model. In fact, 
the latter attempts to take the qualities of each of the former, computing 
homological information in the integer domain avoinding the computation of 
SNF. 

An algebraic minimal model (AM-model) 0] is a tool for computing topo- 
logical information of chain complexes, tool that is applicable to extract 
topological information from digital pictures. In particular, it provides inte- 
ger (co) homology generators, representative (co) cycles of these generators as 
well as the cohomological invariant HBl (derived from the rank of the co- 
homology ring) 0, 0, S]. More concretely, an AM-model for a chain complex 
C = {C, d) is a set ((C, d), (M, d')f, g, 0), where 

(i) M = {Mq} generates a chain complex, denoted by A4, with differential 
d' such that any non-null entry of SNF of the matrix of d' in each dim 
q, is greater than 1. 

(ii) (/, g, 0) defines a chain contraction of C to A^. 

Property (ii) follows that H{C; Z) and H{Ai; Z) have isomorphic homol- 
ogy groups. Moreover, integer (co)homology generators and representative 
(co)cycles of (co)homology generators of C can directly be obtained from Ai. 
See Fig. HJ for example, where Mi = {«,/?}, a being a representative cycle 
of a generator of the free part of Hi{C{S); Z) and /3 a representative cycle of 
a generator of the torsion part. An algorithm for computing AM-models is 
given in 0, S]. It needs to reduce the matrices of the differential to its SNF. 

An algebraic topological model (AT-model) js], [oj is a tool for computing 
over a field Z/p, p being a prime, (co)homology, representative (co)cycles 
of (co)homology generators and the cup product on cohomology of nD dig- 
ital images (and hence also the cohomological invariant HBl). The algo- 
rithm presented in js], jof for computing AT-models over fields runs in time at 
most 0{m^), where m is the number of generators of the given chain com- 
plex. More concretely, an AT-model for a chain complex C = {C,d) is a set 
((C, d),H, f,g,(l)), where 

(i) H = {Hq} generates a chain complex, denoted by "H, with null differ- 
ential. 
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Ho(e(S): Z):^Z 




Figure 4: Homological information obtained from an AM-model for the Klein bottle. 

(ii) {f,g,(j)) defines a chain contraction of C to "H over Z/p. 

Property (ii) follows that the homology groups of C and "H, with coefficients 
over Z/p, are isomorphic. 

Example 3.1. Let K he the simplicial complex shown in Fig. The set 
{{K,d), i^K, /a'5 5'k 5 0k) is an AT-model for K over Z/2, where 

K = {Vi,V2, V3, W4, V1V2, VsVi, V2V3, ViV2Vi}, Hk = {Vi, V2V3} 

and the images of the maps fj^, g^ and (j)^ are: 
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V2V3 
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VIV4 
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VIV2V4 















where a = V2V3 + V2V4 + V3V4. Then, the homology groups of C{K) 
over the field Z/2 are Hq{C{K)]Z/2) = Z/2, Hi{C{K);Z/2) = Z/2 and 
H2{C{K);Z/2) = 0. The representative cycles of homology generators in 
dim and 1 are V4 and a, respectively. 



The main problem of this tool is that from the homology obtained from 
an AT-model, it is not possible, in general, to get all the torsion information 
of the integer homology of the object under study (see Example 12.31) . 
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Figure 5: On the left, the simplicial complex K; in red, the representative cycles of 
homology generators of H{C{K); Z/2). On the right, the elements of Hj^ . 

4. The Notion of A-AT-Model 

As far as we know all the algorithms in the literature for computing 
integer homology (with or without torsion) need to reduce the matrix of the 
differential to its SNF. Our aim is to compute integer homological information 
avoiding the computation of SNF. For this task, we first define the notion of A- 
AT-model which can be seen as a generalization of the concept of AT-model, 
working in the integer domain. We study its properties and give an algorithm 
for computing it. Moreover, we will see that from a A- AT-model for a chain 
complex C, we can obtain the Betti numbers Pg, the prime numbers p involved 
in the invariant factors of the torsion subgroup of H{C;Z), the amount of 
invariant factors that are a power of p, and a set of representative cycles of 
generators of homology mod p, for each p. That is, from a A- AT-model for a 
chain complex C, we can obtain a set X of independent non-boundary cycles 
such that given any cycle c of C, Ac can be expressed as a linear combination 
over Z of the cycles of X. Finally, we present an algorithm for extracting 
this homological information, performed in 0{m^tp{X)) in the worst case, ijj 
being the Euler function. 

First of all, let us define the concept of A-chain contraction as a general- 
ization of chain contraction and study its properties. Later, we will define a 
A-AT-model as a A-chain contraction of a chain complex C to another chain 
complex with null differential. 

Definition 4.1. Let A be an integer, X ^ 0. A X-chain contraction of a 
chain complex C, to another chain complex C with differentials d and d' , 
respectively, is a set (/, g, </>, A) such that: f : C ^ C and g : C ^ C are 
chain maps; fg is Xidc' and (j) : C ^ C is a chain homotopy of Xidc' to gf , 
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that is 



fg = Xidc' 



and 



Xidc — gf = <pd + d(j) . 



Observe that a A-chain contraction with A = 1 is a chain contraction. 
Moreover, if (/, g, 0, A) is a A-chain contraction of C to C with A = — 1 then 
(— /, g, —0) is a chain contraction of C to C. 

Lemma 4.1. // {{C,d), H, f, g,(j), X) is a X-AT-model, with A < 0, then 
((C, d), H, — /, g, —(f), —A) is a X'-AT-model, with X' := —A > 0. 

If {f,g,(j),X) is a A-chain contraction of C to C and (/',(?', 0', A') is a 
A'-chain contraction of C to C", then 



is a AA'-chain contraction of C to C". Another important property is the 
following. 

Proposition 4.2. // there exists a X-chain contraction (/, g, 0, A) (with X ^ 
0) of a chain complex C to another chain complex C then the free subgroups of 
H{C] Z) and H{C'; Z) are isomorphic. As a consequence, the Betti numbers 
of C and C coincide. 

Proof. Let d and d' be the differentials of C and C, respectively. Let a be a 
representative cycle of a homology generator of the free subgroup of H{C\ Z), 
that is, d{a) = and for each /i G Z, /i 7^ 0, yua ^ Imd. Then /(a) is a cycle 
of C since / is a chain map and therefore d'f{a) = fd{a) = 0. Suppose that 
there exists an integer /x 7^ and a chain b' G C such that fif{a) = d'{b'). 
Then, Xa — gf{a) = d(j){a) since {f,g,(f),X) is a A-chain contraction and a 
is a cycle. Therefore, Xfj,a = d{g{b') + fi4>{a)) which is a contradiction since 
a is a representative cycle of a homology generator of the free subgroup of 
H{C; Z). Therefore /(a) is a representative cycle of a homology generator of 
the free subgroup of H{C'; Z). 

Now, let a and b be representative cycles of two different homology gen- 
erators of the free subgroup of H{C] Z). Suppose that /(a) and f{b) are 
representative cycles of the same homology generator of the free subgroup of 
H{C'; Z). Then there exists a chain b' G C such that /(a) - f{b) = d'{b'). 
Then A(a — b) = d{g{b') + (f){a — b)) which is a contradiction. Therefore, /(a) 
and f{b) are representative cycles of two different homology generators of the 



(/7, gg\ A'0 + #7, AA') 



(1) 
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free subgroup of H{C'; Z). Analogously, it is true that if a' and b' are rep- 
resentative cycles of two different homology generators of the free subgroup 
of H{C'] Z) then g{a') and g{b') are representative cycles of two different ho- 
mology generators of the free subgroup of H{C; Z). We can conclude that 
the free subgroup of H{C] Z) and H{C'\ Z) are isomorphic. □ 
Now, define a A-AT-model for a chain complex C as a A-chain contraction 
of C to a chain complex with null differential. 

Definition 4.2. Let C = {C,d)) be a chain complex. Let X be a non-null 
integer. A \-AT-model for C is a set {{C,d), H, f, g, (f), \) , where 

(i) H = {Hq} generates a chain complex, denoted by Ti, with null differ- 
ential. 

(a) {f,g,(f),X) defines a X-chain contraction of C tol-i. 

As a consequence of Prop. 14.21 and since the differential of "H is null, if 
there exists a A-AT-model ((C, d), H, /, g, 0, A) for C then the free subgroup 
of H{C; Z) is isomorphic to "H. 

Proposition 4.3. Given a X-AT-model, a rational AT-model (i.e., an AT- 

model over Q) can directly be obtained as well as rational (co)homology and 
representative (co) cycles of (co)homology generators. Concretely, if{{C, d), H, 
/, g, (j), A) is a X-AT-model for a chain complex C, then ((C, d), H, jf, g, 
is an AT-model for C over Q and {g{h) : h G H} is a set of representative 
cycles of generators of H{C; Q). 



Universal Coefficient Theorem for Homology [IJ, pp. 332] states that 
each Betti number /3g coincides with the vector space dimension of Hq{C; Q). 
In addition, working in the integer domain, a A-AT-model for C provides a 
set of /3q independent non-boundary g-cycles of C over Z, for each q. 

Corollary 4.4. Let {{C,d), H, f, g,(j), X) be a X-AT-model for a given chain 
complex C. Then % is isomorphic to the free subgroup of H{C; Z). Moreover, 
the set {g{h) : h G H} is a set of independent non-boundary cycles ofC over 
Z that generate the set {Xa : a G F{C; Z)}. 

Proof. Let ((C, d), H, f, g, (j), A) be a A-AT-model for C and {hi, . . . , km} a 
set of generators of H. Suppose that {g{hi), . . . , g{hm)} are not independent. 
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Then there exists i, 1 < i < m, such that g{hi) = ^Jli Cjg{hj) for some Cj e 
Z, 1 < j 7^ i < m, and q = 0. Therefore fg{hi) = Yl^=i^jf9i^j)^ ^o that 
= l^jLi CjA/ij. Simphfying, hi = Yl'^=i^j^j^ which is a contradiction. 
Let a be a homology class of the free part of H{C\ Z) and c a representative 
cycle of a, that is, a = [c]. Since /(c) G if, then /(c) = Yl^i^i^iy ^ ^■ 
Therefore, gf{c) = YlT=i(^i9{hi)- So, YILi^^idiK) = Ac - 0rf(c) - #(c). 
Since c is a cycle, then Yl^i^idi^i) = Ac — d(f){c). Considering homology 
classes, Xll^i = A[c]. Consequently, {(yf(/ii), . . . , (7(/i.„)} is a set of 

independent non-boundary representative cycles of generators that generate 
the set {Xa : ae F{C; Z)}. □ 
Similar to rational AT-models, an AT-model over Z/p (where p is a prime 
which does not divide A) can always be defined from a A- AT-model. 

Proposition 4.5. Given a X-AT-model ((C, d), H, /, g, 0, A) and a prime p 
such that p does not divide X, then {{C,dz/p),H, fz/p,gz/p,4>z/p), where 

dz/p = d mod p, fz/p = X~^f mod p, gz/p = g niod p, (pz/p = X'^^cj) f^od p, 

is an AT-model for C over Z/p. Moreover, {gz/p{h) : h G H} is a set of 
representative cycles of generators of H{C; Z/p). 

An important property of A-AT-models is that we can obtain the prime 
numbers p involved in the torsion subgroup of the homology of C. 

Proposition 4.6. Let {{C, d), H, /, g, (p, A) be a X-AT-model. Let a EC such 
that d{a) =0. // there exists h E C such that d{b) = fia where /i G Z, /i 7^ 
and for each p, where < p < fi, pa ^ Im d, then f{a) = and p divides X. 

Proof. Suppose that h E C such that d{h) = pa where p E Z, p and 
for each p, where < p < p, pa ^ Im d. First, let us prove that /(a) = 0. 
Suppose that /(a) 7^ 0. Then pf{a) 7^ (since the ground ring is Z). On the 
other hand, pf{a) = f{pa) = f{d{h)) = 0, since / is a chain map and the 
differential of "H is null. We then have that /(a) = which is a contradiction. 
Now, suppose that p does not divide A, then there exist s, r G Z, such that 
< r < p and X = sp -\- r. On one hand, ra ^ Im d. On the other hand, 
ra = (A — sp)a = Aa — spa = gf{a) -\- (f)d{a) + d(f){a) — sd{b) = d{((){a) — sb) E 
Im d which is a contradiction. We conclude that p divides A. □ 

Corollary 4.7. Let {{C,d), H, f, g, (p, X) be a X-AT-model. If p^(i^p) is an 

invariant factor of Hg{C; Z), then p divides X. 
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Observe that if A is 'large' then there are lots of candidates to take part 
in the torsion subgroup of H{C; Z) (all the primes that divide A). Therefore, 
an important task could be to get the minimum A > such that there exists 
a A-AT-model for C. 

Proposition 4.8. Let C = {C, d) he a chain complex such that in each dim q 
the matrix of dq coincides with its SNF. Let p be the lowest common multiple 
of all the non null elements that appear in the matrices of dq for all q. Define 
the homomorphisms f , g and and the set H as follows: 

• If X E C such that d{x) = fiy for some integer /i > then 

f{x) = , f{y) = , (f){x) = and 4>{y) = - x. 

A* 

• IfxEC such that d{x) = and there is not any z E C being d{z) = fix 
for some integer /i > 0, then 

f{x) = X , g{x) = px and (j){x) = . 

• H = {x E C such that d{x) = and there is not any z E C such that 
d{z) = px for some integer p > 0}. 

The set {{C,d), H, f, g,(j), p) defines a p-AT-model. Moreover, this p-AT- 
model satisfies that a prime p divides p if and only if it takes part in the 
torsion subgroup of H{C] Z). 

Observe that in order to obtain a A-AT-model for a chain complex C with 
minimum A, we need to compute the SNF of the matrix of the differential of 
C in each dimension. The following algorithm computes a A-AT-model for a 
chain complex C without computing the SNF of the differential. In this case, 
A might be non-minimum. 

Algorithm 4.9. Computing a X-AT-model for a chain complex C. 
Input: A chain complex C = {C,d) of dim n. 

Hq:=Co, a := Co , A := 1 , /o := idc^ , go ■= idno > 0o := . 
For q = 1 to q = n do 

Hq:={}, x:=0, ^ := , ^ := . 
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While M := {\{fg^idq{c),c')\ 7^ : c e Cg, c' e i^g-i} is not empty do 
Take a G Hg^i, f3 e Cg s.t. | (/g„irfg(/3), a) | =min M then 
Hg^,:=Hg^,\{a}, A:=AU{P}, /,(/?) :=0, := 0, 

X := a), 7 := 7' := A/3 - . 

For each b E Cr , 0<r<g — 1, do 

Mb) ■.= xfr{b)-{fr{b),a)^, 

(f)r{b) ■.= x(j)r{b) + {fr{b),a)i , 
\ := x\. 
For each a E Cg \ A do 

Hg := Hg U {a}, A := A U {a} , 

/g(a) := a, gg{a) := Aa - (j)g^idg{a), (j)g{a) := 0. 

Output: The set {{C,d),H, f,g,(j),\) . 

Fixed q, there are two different parts in the algorithm. First, while the 
matrix M corresponding to fq-idg is non-null, take a G Hg_i and (5 E Cg 
such that I a) I is the minimum of M (without considering signs). 

Then, a is deleted from Hg_i (i.e., a homology class is destroyed) and /3 is 
added to A (i.e, /3 has been used). Second, if the new redefined matrix M is 
null, add, to Hg, the non-used elements of Cg (i.e. new classes of homology 
are created). At the end of the algorithm, the set A consists in an ordered 
set of all the elements of C. 

Theorem 4.10. The set ((C, rf), if, /, (7, 0, A) obtained applying Alg. |^.^ de- 
fines a X-AT-model for the chain complex C = {C,d). 

Proof. Let A = {a^, . . . , a*"} such that if z < j then a* was added to A before 
when the algorithm was applied. Fixed i, assume that fr-idr{a^) = 0, 
Xa'^ — grfrioJ') = (t>r-idr{a'') + dr+i(j)r{a'^), fr+i4>r{<^^) = and, if G Hj., 
drgr{o!') = and frgr{o!^) = Xa^, ioi 1 < k < i and r = dim{a^). 

Let q = dim{a^). We have to prove that the redefined maps that we will 
denote here by /', g' and 0' and the new integer A' satisfy that fl._idr{a^) = 0, 
X'a'' - g'^fl.{a^) = (j)'^_^dr{a^) + rfr+i0^(a''), /^+i0^(a'') = and, if G H^, 
drg'ria'^) = and flg'j.{a}^) = X'a^ for 1 < < i and r = dim{a^). 

We have to consider two cases. (1) If a* was added to A while M was 
non-empty, then g' = g (that is, g does not change). Therefore, if G Hr 
for some r, then drg'ria^) = drgria^) = for all k, 1 < k < i. Consider the 
element a G Hg_i that defines the auxiliary parameter x = {fq~idg{a^) , a) at 
this stage. 
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If k = i, then a* ^ Hg. We have that fg_idq{a'^) = xfq^idq{a^) — 
{fg^idq{a'),a)'y = x'j - x-f = and 0g_i(i<^(a*) + rf,+i0^(a*) = ^^.^^^(a*) = 
x4>q-idq{a') + (/g_idg(a*), a)7' = x{(j)q-idq{a') + Aa* - (f)q^idq{a')) = xXa' = 
A'a*. Moreover, /g+i0g(a*) = since 0^(a*) = 0. 

Observe that: For any r-chain b such that r < q — 1, f^{b) = xfr{b) 
and 4>r{b) = xcprip). dq{^') = dq{Xa^ — 0q_ic/g(a*)) = Xdqi^a") — {\dq{a'^) — 
gq-ifq-idq{a') - (t)q_idq_idq{a')) = gq-ifq-idq{a') = Qq-iii). For any g-chain 
6, fq{b) = and (f>q{b) = 0, therefore, fq{b) = and 0^(6) = 0, since a* was 
added to A while M was non-empty. 

For !</;;<«. (l.a) If r := dim{a'') < q — 1 then fl_idr{a!^) = 
xfr~idr{a^) = by induction and (f)[._idr{a'') + dr+i(pr{a'') = x{(j)r~idr{a'') + 
dr+iMa'')) = x{\a'' - grfr{a^)) = \'a^ - grfM^). Moreover, /;,+i</);(a*^) = 

X fr+l4>r{.(i'') = x{xfr+l(pr{(l'') — {fr+l4>r{(l''),C()l) = by iuductioU. If fl'^ G Hr 

then flgria}') = xfr-gria^) = x\a^ = X'a''. (l.b) If di7n{a^) = g — 1, 
then fq_2dq_i{a^) = xfq^2dq-i{a^) = by induction and 4>'q^2^q-ii^'') + 

dq(j)'q_J^a^) = X(t)q-2dq-l{a^) +dq{x(f)q-l{a'') + {fq-l{a'') , a)'y') = X{(f)q-2dq-l{a'') + 

dq(j)q-i{a^)) + {fq-i{a''), a)dq{i) = x{\a^-gq-ifq-i{a^)) + {fq-i{a!'),a)gq-i{-i) = 
\'a^-gq.,rq_,{a^). Moreover, fqCj^'q.M) = rq{xct>q.,{ a^) + {fq^M).Oi)i). 
Now, f'q(t)q_i{a^) = since for any g-chain b, f^{b) = 0; fgi-y') = f^{Xa' - 
(t>q-idq{a^)) = for the same reason than above. Then, f'q4>'q_i{a^) = 0. If 
a'' e Hq^i then a'' ^ a and fq_igq^i{a^) = xfq^igq^i{a'')-{ fq^igq^i{a^), a)7 = 
xXa!' — (Aa'^,a)7 = x\a^ = X'a!'. (l.c) If dim^a!') = q, then fg_^dq{a!') = 
X fq-idq{a'') — {fq^idq{a^) , 0)7 = by induction and (j)'q_idq{a^)+dq+i(l)'q{a^) = 
x{(l)q_idq{o!') + dq+i(t)q{a!')) + {fq-idq{a'') , a)-f' = x(Aa^ - gqfq{a^)) = X'a'' - 
gqfg{a^). Moreover, /^+i(/>^(a'=) = 0. 

(2) If a* was added to A when M was empty, then A' = A (that is, A 
does not change). First, f'q_idq{a^) = fq-idq{a^) = by induction and Aa* — 
= Xa'-g'g{a') = Aa*-(Aa*-0g_irfq(a')) = (j)q-idq{a') = (f)'g_^dq{a') = 
4>'q-idq{a'') +dq+i(f)'g{a^) . Second, dqg'g{a^) = Xdq{a^) — dq(f>q^idq{a^) = Xdq{a^) — 
iXdq{a')-gq_i fq_idq{a')-(f)g^dq_idq{a')) = gq-^ifq-^idg{a') = and fl^g'g{a') = 
fg{Xa^ — (t>q-idq{a^)) = Aa* — fq(pq-idq{a^) = Xa\ since fr+i(pr{ci^) = for 
1 < k < i and r = dim[a^), by induction. Finally, = since 

0;(aO = 0. □ 

To study the complexity, consider the order of the elements of C given 
by the set A. Fixed z, count the number of elementary operations involved. 
We have to update frip) and (prip) for each 6GCr,0<r<g— 1. Observe 
that in order to update fr{b) and 0r-(^) in the worst case, we have to add to 
both a chain containing generators (if is the number of generators of 
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Figure 6: A visual interpretation of the maps fs and 0s on a triangulation of the Klein 
bottle when applying Alg. 14.91 a) after adding all the vertices and edges; b) after adding 
all the vertices edges and triangles except for v^vqVs] c) at the end. 



dim q). In the worst case, rriq could be of the same order than m (the total 
number of generators of C). So, fixed i, the total cost of these operations is 
0{im?). Therefore, the total algorithm runs in time at most 0{m?). 

Example 4.11. Consider the simplicial complex S derived from the trian- 
gulation of the Klein bottle given in Fig. [7] and the chain complex C{S) 
of dim 2 associated to S. Applying Alg. \4.9[ we obtain a 2-AT-model, 
{C{S), Hs, fs, gs, 0s, 2), for C{S). For the sake of simplicity in this example, 



each time a negative value for A appears, we apply Lemma 4jJ_ one? turn it 
into positive. After running Alg. \4-S\ only for the vertices and edges of S, we 
get that, at this stage, = {vq} U {x : x is an edge marked in red in Fig. 
^a)} and the values of fs, gs o-nd 0s at this stage are: 





Vo 


Vi 
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X 


fs 


Vo 


Vo 





X 


9s 


Vo 






X - (ps'diix) 


4>s 














where I < i < 8, y denotes an edge drawn in black in Fig. \^a) and 'y{vo,vi) 
is a path connecting the vertices vq and Vi, i ^ following the arrows drawn 
in black in Fig. \^a); for example, ■y{vo,v6) = VqVi + V1V2 + V2Vq. 

Now, apply the algorithm for all the vertices, edges and triangles except 
for the triangle v^v^vg ( to get a better idea of the geometric meaning of the 
maps fs, 9s o,nd 0s, suppose that we do not consider this triangle until the 
end). We get that Hg = {i^o, ^0^3}- The values of fs, 9s and 0s at this 
stage are: 
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where z denotes one red edge (see Fig. ^h)) andt a triangle except for v^v^v^. 
7^ is 0, -^VqV^, -^VqV2 or VqV2 — v^v^ if the arrow drawn over z in Fig. ^h) 
is, respectively, red, blue, green or black. 7^'^-) is a 'path' from z following the 
arrows drawn in Fig. \^b). For example, 'y(vev7) = v^vqVj — v-^v^Vq + ViV2V^ + 

ViVr^Vj - - V0V3V5 + V0V1V5 + VoV^Vy - VqViVt. 

Now, add the triangle V4^vqVs. Then fsdiv^v^v^) is ^vqv^. Therefore, the 
new images for fs, Qs (^nd 0s ore; 
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V0V2 






t 




fs 


2vo 







2VoV2 













9s 


vo 






V0V2 - V1V2 - VqVi 










4>s 
























where 7^," is or ±VoV2 if the arrow drawn over z in Fig. \^c) is, respectively, 
red or green. 7^^^ is a 'path' from z following the arrows drawn in Fig. \^c). 



The following proposition shows that AT-models over Z/p, where p is 
prime, can also be computed using Alg. 14.91 

Proposition 4.12. Working with coefficients in Z/p, where p is prime, let 

(C, dz/p), Hz/p, fz/p, dz/p, 4>z/p, A) 

be the output of Alg. \4.9[ Then 

{{C, dz/p), Hz/p, X ^fz/p,g,X Vz/p) 

is an AT-model over Z/p. Furthermore, {gz/p{h) : h G Hz/p} is a set of 
representative cycles of generators of H{C; Z/p). 

In order to obtain a A-AT-model with A 'small', the following algorithm 
can be used as a preprocessing of Alg. 14.91 It only considers the case in 
which {d{a'),a) = 1. The output of this algorithm is a chain contraction of 
the given chain complex C to another chain complex C with less number of 
generators. We can compute a A-AT-model for C applying Alg. 14.91 to C. 
We can obtain a A- AT model for C, composing the chain contraction of C to 
C with the A-AT-model for C obtained before, using Formula ([1]) given in 
pagefTTl 

Algorithm 4.13. Preprocessing. 
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Input: A chain complex C = {C,d) of dim n. 

Cq := Co , d'Q = do, A := Cq , fo := idc^ , go ■= idc^, , <Po ■= ■ 
For q = 1 to q = n do 

C'^:={], x:=0, 7:=0, i := Q . 

While mm{\{fq-idq{c),d) \ ^ : ceCq.c'e C^_J = 1 do 

Take a G q„i, 13 eC^ s.t. | a) | = 1 then 

:= \ {a}, A := A U {/?}, /,(/?) := 0, := 0, 

a; := {fq~idq{l3), a), 7 := fq-idq{l3), y ■= /3 - (f)q_idq{l3) . 
For each b G Cg_i do 

:= fq~i{b) -x(/g_i(6),a)7, 
:= +x(/g_i(6),a)7', 
For each a G \ A do 

C^:=C;u{a}, A:=AU{a}, 

o?^(a) := fq-idq{a), fq{a) := a, ^(^(a) := a-0g_idg(a), 0g(a) := 0. 
Output: The set ((C, rf), (C'.ti'), /,^, 0, A) . 

Theorem 4.14. The set ((C, d), (C, d'), /, ^f, 0) obtained applying Alg. \4.13 
defines a chain contraction of the chain complex C = (C, d) to the chain 
complex C = (C, d'). 

The proof of this theorem is given in the appendix. 

5. Extracting Integer Homological Information 

As showed before, a A-AT-model for a given chain complex C provides 
information of the free subgroup of if (C; Z) as well as the prime numbers p 
involved in the invariant factors of H{C; Z). In order to obtain representative 
cycles mod p of generators of the free and the torsion subgroups of H{C] Z) 
we have to compute AT-models for C over Z/p, for each prime p dividing A. 

Algorithm 5.1. Computing integer homology information and representa- 
tive cycles of homology generators of a chain complex C. 

Input: a chain complex C = {C,d) of dim n. 
Apply Alg. 14 . 91 with coefficients in Z for computing a 
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A-AT-model {{C,d), HJ, g, (p, \) for C; 
For q = to q = n do 

Pq := number of elements of Hg] 
G:={g{h): h e H} . 
For each prime p dividing A do 

Apply Alg. 14.91 with coefficients in Z/p to compute an 

AT-model for C over Z/p: ((C, dz/p), ^^z/p, /z/p, fi-z/p, 0z/p) ; 

Pii^p) '■~ number of elements of ifz/p in dim q. 

T{o,p) = /3(o,p) - ; 

For g = 1 to n do 

Cz/p = {9z/p{h) : h G ifz/p} • 
Output: The sets G, . . . ,/3„,}, 

{Gz/p : p is a prime dividing A}, 

and p) : < q < n and p is a prime dividing A} . 

After computing Alg. 15. II for a given chain complex C, we obtain: 

• the Betti numbers /3g for < q < n, and a set G of independent non- 
boundary cycles of C over Z (in fact, G is also a set of generators of 
HiC;Q)y, 

• the prime numbers p involved in the invariant factors corresponding 
to the torsion subgroup of H{C, Z), the amount of invariant factors in 
each dim q that are a power of p, T(^q,p), and a set G^/p of representative 
cycles of generators of H{C;Z/p), for each prime p dividing A. 

This algorithm uses Alg. 14.91 over Z for computing a A-AT-model for C 
and over Z/p for computing AT-models for C, for each prime p dividing A. 
Since the complexity of Alg. 14.91 is O(m^), then the complexity of Alg. 15.11 
is 0{m^ip{X)) in the worst case, ip being the Euler function. 



Example 5.2. In Example 4-11. we have applied Alg. \4.9 and computed a 



2- AT-model for C{S). We obtained that Hg = {vo,vqV2} and then the Betti 



numbers of C are /3o = 1, /3i = 1 and (32 = 0. Now, apply Alg. with 
coefficients in Z/2 to obtain an AT-model ((C, (iz/a)), /z/25 fl'z/25 0z/2) 
where Hz/2 = {t'o, ^o^'2, ^^0^4, t'o^'s^'s} forC{S). Then, /3(o,2) = 1, /3(i,2) = 2, 
/3(2,2) = 1- Therefore, T(o,2) = 0, T(i,2) = 1 and T(2,2) = 0. 

We conclude that Ho{S; Z) ~ Z and Hi{S; Z) ~ Z © Z/2. 



20 



Figure 7: A porous 3D digital image and representative cycles of its 2709 black holes. 



Alg. IS.ll has been implemented |9|. Fig. [5] shows a 3D digital image with 
18491 black voxels. In 144 seconds we obtain that the image has 51 black 
connected components, 2709 black holes and 8 black cavities. 

6. Conclusions and Future Work 

We have constructed an algebraic tool, called A-AT-model, in order to 
compute topological features of n-dimensional objects (for instance, nD dig- 
ital pictures). Associating an algebraic representation to the initial object, 
a chain complex, the algorithm presented here provides the free subgroup 
of homology, with integer coefficients, as well as the primes p that are the 
candidates to be involved in the torsion subgroup. In fact, such a set of 
candidates is the one of divisors of A. We have determined the minimum A 
valid for the existence of a A-AT-model. We have also presented an algorithm 
which gives a A-AT-model with 'smaller' A than the one given in [loj. As for 
the complexity, the algorithm presented is performed in 0(m^'i/'(A)) in the 
worst case, ip being the Euler function and m being the number of cells of 
a cell decomposition of the object. In order to improve the complexity, first 
compute a chain contraction from the initial chain complex using Alg. 14.131 
to a 'smaller' one with same integer homology and then compute the A-AT- 
model for the latter. Afterwards, in order to obtain the amount of invariant 
factors that are a power of p and representative cycles mod p of homology 
generators, for each p we compute an AT-model with coefficients in Z/p. 

Future work can be addressed to obtain generators of the torsion sub- 
group of H{C;Z), with integer coefficients. Another task is to study the 
computation of cohomology features over Z from a A-AT-model. 
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Appendix: Proof of Theorem 14.141 

Let A = {a^, . . . , a™} such that if z < j then a* was added to A before 
when the algorithm was applied. Fixed i, assume that fr-idr{a'') = d'j.fr{a''), 

— grfrio}') = 4>r-ldr{ci^) + C^r+lf^r (o'^) , fr+l4>r{ci^) = and, if tt^ G C^, 

drQriO''') = gr-id'rQriO'^) a^d frdriof^) = , foT 1 < k < 1 and r = dim{a^). 
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Let q = dim{a^). We have to prove that the redefined maps that we 
will denote here by d" (the redefined map of c?'), /', g' and 0' satisfy that 

and, if G C',, drg'^o^^) = g'^_id'^g'^{a^) and /^fi'',(a'') = a'', for 1 < A; < i 
and r = dim{a^). 

We have to consider two cases. (1) If a* was added to A while mm{| (/q_i(iq(c), 
c')| 7^ : c e Cq, c' e Cg„J = 1, then d" = d' and g' = g (that is, d' and 
do not change). Therefore, if G Hr for some r, then drg'j.{a^) = drgria^) = 
gr-Ma^) = g'_ id'liO'k) for all 1 < k < i. Consider the element a e C'q_i 
that defines the auxiliary parameter x = {fq-idg^a^) , a) at this stage (observe 
that X = ±1, then xx = 1). 

lik = i, then a* ^ C'^. We have that fq_idq{a^) = fq^idq{a^)—x{fq^idi{a'^), 
a)7 = 7 - = = dqfq{a') and ^'^.^^^(a*) + 0^5+10^(0*) = (j)'q_^dq{a') = 
(l)q-idq{a^) + x{fq-idq{a^),a)i = (pq-idq{a')+xx{a'-(j)q-idq{a'-)) = a\ More- 
over, fq+^^'q{a') = since ^'q{a') = 0. 

Observe that: For any r-chain b such that r < q — 1, f',.{b) = fr{b) and 
0^(6) = (j)r{b). dq{i) = dq{a'-(j)q.idq{a')) = dq{a')-{dq{a')-gq^ifq.idq{a')- 
(j)q-idq_idq{a')) = gq-if q-idq{a') = gq-ii'j). FoT auj g-chain b, fq{b) = and 
4)q{b) = 0, therefore, /^(6) = and 0^(6) = 0. 

For 1 < k < i. (l.a) If r := dim{a'') < q — 1 then f^_^dr{a'') = 
fr-idr{a^) = d'^_^fr{a^) = d'^__^fl{a^) and (t)'^_idr{a^)+dr+i(j)'^{a^) = (j)r-idr{a'') + 
dr+iMa'') = a^-grfr{a^) = al'-grfria!")- Moreover, f^^M(^'') = fr+iMa'') = 
/r+i0r( o.^) —x{fr+i4>r{ci^), Ci)l) = by iuductiou. If G C', then flgr^o!^) = 
frgAaf') = a^- (1-b) If dim{a'') = g - 1, then f'q_^dq-i{a}') = fq-2dq-i{al') = 

d'q_Jq^l{a^) = d'q_Jq^l{a'') - x{f q^l{b) , O) f q-idq-ldq^^) = d'q_Jq-l{a'') - 
d'q_,{x{fq^l{b), a)fq^2dq{f3)) = d'q_J'q_,{a^) = d^.J'q.^) ^nd ct>'q^,dq^r{a^) + 

dq(t)'q_i{a^) = <Pq-2dq-i{a^)+dq<pq_i{a^)+x{fq.i{a^),a)dq{i) = a^-gq^ifq^i{a^) + 
x{fq_i{a^),a)gq_i{-i) = a^-gq-ifq_i{a^). Moreover, f'q(t)'q_^{a^) = /^(0g_i(a'=) + 
x{fq-i{a^),a)'-)'). Now, /g0q_i(a'^) = since for any g-chain 6, f'^ib) = 0; 
fqil') = fqi^^ ~ <Pq-idq{a^)) = for the same reason than above. Then, 
/^0^_i(a'=) = 0. If a^' G C^„i then a^' ^ a and /^_i^g„i(a^) = /g_i^g_i(a'=) - 
x{fq^igq^i{a^),a)j = — x{a^ ,a)'~^ = a^. (l.c) If dim{a^) = q, then 
fq_idq{a^) = fq-idqial') - X {fq_idq{a^) , a)^ = d'Jq{a^) - x{d'Jq{a^),a)-f = 
= d'l^f'q{a^) and (t)'q_M(^^) + = 0Q-ic?g(a'') + dq+i(t)q{a^) + 

x(^_icig(a'^),a)7' = a'^ - ^Jg(a'^) + x{d'qfq{a^),a)i = a'' - gqfq{a^) = 
- dqfqia'')- Moreover, f'q+i<p'q{a^) = 0. 

(2) If a' was added to A when {\{fq_idq{c),c')\ ^ : c e Cq, c' e 
C^„J = or rmn{\{fq^,dq{c),c')\ ^ : c E Cq, c' E C'^_,} > 1 then, 
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= fg-id,{a^) = <(aO = </^(aO and a\- g'J^{a^) = a' - 
g'^(a') = a' - (a* - 0g_ic/g(a*)) = 05_irfg(a*) = ^'^.^^^(a*) + rfg+i0'g(a'). Sec- 
ond, dqgg{a') = dq{a') - dq(j)q-ldg{a') = gq-lfq-ldq{a') + (f)q-2dq-ldq{a') = 

gq^Jq^,dq{a^) = ^?,-i<(aO = g'q.d'^ia') and /^(aO = fqia^ - ^,-idq{a')) = 
~ fq(pq-idq{a') = a\ since fr+i4>r{,o!^) = for 1 < A; < z and r = dim{a^), 
by induction. Finally, fq^i4>q{a'') = since = 0. 
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STATEMENT 



We present some notes in response to the suggestions made by the re- 
viewers. 

Concerning the comments of Reviewer 1, the main remark is that we 
should add further explanations for the readers who are not specialist in the 
field. We have followed this philosophy. We could classify the suggestions in 
the following categories: 

• Little gaps or changes of notations: we have followed all the reviewer's 
instructions. 

• Addition of examples: we have added more examples to illustrate some 
aspects that the reviewer asked to be clarified in this way. 

• Further explanations: we have recalled more background concepts and 
detailed some technical aspects at the points that the referee suggested. 
In some other cases, we have simply extended the given explanations. 

• Regarding Algorithm 1 (Algorithm 13 in the previous version), it has 
been modified since it was difficult to follow the technical aspects, as 
pointed out by the reviewer. The changes have allowed to carry out the 
proof of Theorem 16 (Theorem 14 in the previous version) in a clearer 
way. 

• Proof of Theorem 19 (Theorem 18 in the previous version) has been 
provided in an appendix. 

• With respect to the reviewer's question about the reduction of the 
complexity using Alg. 2 (Alg. 17 in the previous version), we have to 
say that it is difficult to quantify such as a reduction in general. We 
plan to study this question in detail and expect to provide experiments 
and concrete examples of such a reduction. 

As for the comments of Reviewer 2, all the suggestions have been carried out. 
Most of them are Grammar errors or little gaps that have been corrected. 
Also the proposed change in the numbering for definitions, theorems, etc., 
has been undertaken. 
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